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ABSTRACT 

We use the phase space position-velocity (x, v) to deal with the statistical prop- 
erties of velocity dependent dynamical systems, like dissipative ones. Within this 
approach, we study the statistical properties of an ensemble of harmonic oscillators 
in a linear weak dissipative media. Using the Debye model of a crystal, we calculate 
at first order in the dissipative parameter the entropy, free energy, internal energy, 
equation of state and specific heat using the classical and quantum approaches, for 
the classical approach we found that the entropy, the equation of state, and the free 
energy depend on the dissipative parameter, but the internal energy and specific heat 
do not depend of it. For the quantum case, we found that all the thermodynamical 
quantities depend on this parameter. 
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1. Introduction 



Our foundation of classical and quantum statistical mechanics [1] is based on the 
Hamiltonian formalism [3] of conservative systems. These systems are particular 
cases of much more general ones called autonomous systems which are those where 
the total force acting on the particle does not depend explicitly on time, other- 
wise they are called non-autonomous. For conservative systems, the Hamiltonian 
is a constant of motion of the system, and, in principle, one can calculate all the 
thermodynamic characteristics of an ensemble of N particles governed by this Hamil- 
tonian, through the acknowledge of its associated partition function. For dynamical 
systems which depend explicitly on the velocity, for example dissipative systems. 
There are mainly two approaches to study these systems. The first one uses the 
kinetic equation (Vlasov or Boltzmann equations) to find the density distribution 
and, then, to calculate the desired statistical property. The other one tries to find 
an effective or phenomenological Hamiltonian and, then, to estimate the desired 
thermodynamic quantities from the stationary statistical mechanics (statistical en- 
semble) [5] . This last approach is possible since for autonomous systems sometimes is 
possible to find time-independent constant of motion, Lagrangian, and Hamiltonian. 
However, it is known that the Lagrangian and Hamiltonian formulation have some 
problems. First, the Lagrangian (therefore the Hamiltonian) may not exist for some 
dynamical systems [6]. Second, There could be two completely different Hamiltoni- 
ans which describe the same classical mechanical system, but completely different 
quantum and statistical systems [7]. Third, time-explicitly dependent systems have 
an ambiguous Lagrangian and Hamiltonian formulation [8]. Fourth, Ambiguity La- 
grangian and Hamiltonian formulation happens even for the harmonic oscillator [9] 
(where a problem of consistency of units arises). Finally, there are autonomous 
systems where one can not have explicitly the Hamiltonian of the system because 
the inverse relation v = v(x,p) can not be given from the original definition of 
p = p(x,v) [10]. This last problem happens particularly on dissipative systems and 
is the main reason one would like to have an alternative approach in mechanics [11], 
statistical and quantum mechanics to deal with them. 

In this paper, we take the original ideas of Maxwell in statistical mechanics and 
Heisenberg in quantum mechanics to study these areas from the point of view of 
coordinates and velocities rather than coordinates and generalized linear momen- 
tum. We apply this approach to the ensemble of linear oscillators within a linear 
dissipative media, and we will do this from the phenomenological point of view and 
at first order in the dissipation parameter in the constant of motion. 

2. Dynamical equation in the space (x,v) 

In this section, we restrict ourselves to one-dimensional single particle motion to 
show the main ideas. The approach will be extended immediately to N independent 
particles moving in a three-dimensional space. The equations of motion of a particle 
moving under a time-independent force can be written as the following autonomous 
dynamical system 

x = v , (la) 
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and 

v = l^l, (lb) 
m 

where m is the mass of the particle, x and v represent its position and velocity, and 
f(x,v) is the total force acting on the particle. A constant of motion of this system 
is a function K(x, v) such that dK/dt = 0, that is, it satisfies the following equation 

ox m ov 

Given this constant of motion, the Lagrangian of the system is given by [15] 

L = v [ ^f^dv + A(x)v , (3) 



where A(x) is an arbitrary function, and the term A(x)v represents the gauge of the 
Lagrangian. Thus, the generalized linear momentum is given by 

P =d^^ (4) 

and it is here where our analysis starts. Suppose now that it is not possible to get 
v = v(x,p) from (4). Then, the Hamiltonian of the system, H = vp — L(x,v), can 
not be given explicitly, and it will be given in implicit form through the constant 
of motion. Therefore, this constant of motion (with units of energy) can help us to 
avoid this an other problems already mentioned in the introduction. 

Consider N particles moving in the three-dimensional space and under a time- 
independent force acting on them (autonomous system). The motion of these 
particles is restricted on the hypersurface of the space !ft 6Ar defined by the time- 
independent constant of motion K(pc, v), where x, v € 5R 3Ar . Thus, for a canonical 
ensemble of iV-particles , the classical partition function would be given by 



Z "' 



h m T] 



J e -/^(x,v) dx(fv ; (5) 



where m is the mass of the particles, h is the Plank's constant, r/ = 1 for distin- 
guishable particles and rj = N\ for non distinguishable particles, f3 = 1/kT with k 
being the Boltzmann's constant and T the temperature of the system, and tixtiv is 
the measure in the space 

3N 

(ixiiv = Y[ dxjdvj . (6) 

3=1 

For a quantum canonical ensemble of iV particles in a batch temperature T, the 
quantum partition function would be 

Z = ]>>e-^, (7) 

i 

where uji represents the degeneration of the eigenvalue where this eigenvalues 
comes from the solution of the equation 

K(St, v)ipi = Eitpi . (8a) 
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K, x and v are the Hermitian operator associated to the constant of motion K(x, v), 
the position x, and the velocity v. ^i(x) is the eigenfunction which is related with 
the wave function ^(x, t) of the Schrodinger equation, 

^^l^ = K(x,v)* i (x,t) , (86) 



as 



*(x,t) = ^Qe-^V i (x) , (8c) 



where \Ci\ 2 represents the probability that the system be on the state tpi. The 
position and velocity operators are defined and related by the following expressions 

Xj = xj , Vj = -— , [xi,Vj] = —5ij, [xuXj] = [vi,Vj] = . (8d) 

m oxj m J 

In this way, with Eq. (5) and Eq. (7) it is possible, in principle, to study the 
thermodynamic characteristics of an ensemble of N particles of any autonomous 
system through the classical or quantum canonical partition functions. Of course, 
this approach is reduced to that one of the Hamiltonian formalism whenever this 
one exists explicitly [11]. Note that the same approach can be made for any other 
type of ensemble. 



3. Constant of motion 

In reference [12], a constant of motion K a (x,v) was given for the dynamical system 

x = v (9a) 

and 



v 



2 a 

-to X V , 

m 



(9b) 



where uj is the free natural frequency of oscillations, m is the mass of the particle, 
and a is the coefficient of the dissipative force which arises phenomenologically as 
an average effect of the interaction with the particles of the medium. The constant 
of motion of this system is given by 



K a [x, v) = — I v + 2u> a xv + u> x \e~ 



■2uj a G(v/x,uj,uj a ) 



where uj a and the function G are defined as 

a 



(10a) 
(106) 



and 



i 



In 



2 - ^ U+£+ v^-^ 2 
i 

l 



: arctan 



uj 2 - UJ 2 a 



if uj 2 < io\ 



if uj 2 = u? a 



if to 2 > uji. 



(10c) 
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This constant of motion has the following limit 

lim K a (x, v) = Kq(x, v) = -mv 2 + -muj 2 x 2 . (11) 
atoO 2 2 

For very weak dissipation (uj a <C uj) and at first order in oj a , the constant of motion 
can be written as 



K(x,v) = K (x,v) + — 

UJ 



mujxv — 2Kq(x, v) arctanl 



\ujx J 



(12) 



In principle, one can get the Lagrangian of the system through the expression (3), 
and from this Lagrangian one could get the generalized linear momentum of the 
system p = p(x,v). However, it is not possible to know v = v(x,p), even for the 
very weak dissipation case. Thus, the Hamiltonian is given implicitly through the 
above constant of motion. 

The quantization of (12) was carried out in reference [13], where the modification of 
the eigenvalues of the harmonic oscillator (En^ = fiuj(n + 1/2)) were given at first 
order in perturbation theory as 



E n {u) = Huj (n + ^ (\ - 



TTUJc 
UJ 

oo 21+1 , i/o\2/oiii\2 



( 2 n + l ) _gg7-V2Vr^ + 1 V (2n-2i- a ) 



^ 3 4 tif^^ 1 J ^ s ' (2/ + l) 2 (2^ + l + s 

(13) 

where the correction has been made up to second order in uj a . 
4. Classical approach 

Consider an ensemble of N particles moving independently which their equation of 
motion is given by the following dynamical system 



where ujj is the natural free frequency of oscillations. The dissipative coefficient 
is the same of all oscillators since they are inside the same medium. One can 
think, for example, in a crystal which is inside a bosonic or fermionic medium 
where the wave length of the particles of this medium is less than the separation 
of the crystal components [16]. The average interaction of the crystal components 
with the particles of the medium may be modeled by Eq. (14). In Astrophysics, 
one could think about the core of a neutron star as a superconducting lattice of 
protons inside a dissipative medium generated by the huge amount of neutrinos 
(antineutrinos) appearing from the weak decay of neutrons and protons. Since the 
density in the core of these stars is very big (10 15 gr/cm 3 )[17] the interaction with 
neutrinos (antineutrinos) is not so neglectful. Thus, maybe the local motion of an 
ensemble of nucleons in the core could be model by Eq. (14). 

Let us see first the implications of using classical canonical partition function to 
study the thermodynamical characteristics of the system. To do this, let us make 
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the following change of variable 



x« = * — — cos d>j , and vi = \ — sinfffl,-) , (15) 

y mcjj y m 

where the variables <f>j and J, have the following variation <fij £ [0, 2ir] and J,- € 
[0, oo). Then, from (12), the jth-constant of motion and measure are given with this 
coordinates as 

Kj = LOjJj — uj a Jj(sm2(pj — 2(f) j) , (15a) 

and 

dxjdvj = — dJjdcpj . (156) 

In this way, the partition function (5) is written as (77 = 1 since the particles in the 
crystal are distinguishable) 

37V 37V . 3JV 1 . 

Z = T^tv II / e-PKste'^dXjdvj = Hl e-Mte-^^Wj-WdJjdtj (16) 
j=i j'=i 

which can be integrated, bringing about the following result 

*=n^. dr.) 

where the function /(£) has been defined as 

Fig. 1 shows the behavior of this function as a function of oj a /uj. Note that 1(0) = 
2tt, and that w Q /w must be less than the unit, according to our approximations of 
expression (5). From (17a), one gets 

37V f j, I 



37V / 

\nZ = -3iVln/i + ^ln 

7=1 V 



Since iV is big and there are 3iV normal modes of oscillations in the crystal, one can 
assume continuity in the frequency spectrum and write (18) as 

In Z = U + I" In (l^±M) g ( u ) di0 , (i 9a ) 



f3uj 

where Uq has been defined as Uq = —3N\nh, and g(u) is the density spectral which 
must satisfied the condition 

roo 

/ g(u) dco = 3N . (196) 
Jo 

Using the Debye's model of solids [14], the spectral density is given by 



ffM = < 



9iYcZ ifO<a;<u, D 



w fl , (20a) 

, if LO > Li>£) 
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where wo is the Debye's frequency of the solid which is defined by the cutoff fre- 
quency such that 



and it is given by 



PLU D 

Jo ' 



UJ D 



g(uj) duo = 3N . 



3N \ V 3 
J v c 



4irV 



(20b) 



(20c) 



The variable V represents the volume of the solid, and v c is the average velocity of 
the elastic waves in the solid which is given in terms of the longitudinal (vi) and 
transversal (vt) waves as 

3 2 1 



v c vf vf 



Substituting (20a) in (19a), it follows that 



9iV p 2l (I(uj a /uj)\ 
In Z = U Q + — J uo 2 In ( v y 1 } duJ 



D 



(20d) 



(21) 



In this way, one can calculate the thermodynamic characteristics of the system. The 
internal energy of the system, its entropy, its specific heat, its equation of state, and 
its free energy are given by 



I -^ = 3NkT, 



S = k 



U + 



9N 



D 



P 



3NkT 
V 



In 



C V 

l(UJg/uJ D ) 

(3uj d 



dp 



uj 2 In 



df 
9N 



I(ulg/Uj) 

pLO 



3Nk 



9N ["d . 



: ln 



duj + 3N 



I{uj a /uj) 



Plo 



doj 



and 



-kT 



rr 9N 

iO 



UJ 



D 



UJ' 



; ln 



I{uJ a /uj) 
Puj 



duo 



(22a) 



(226) 



(22c) 



(22d) 



(22c) 



As one can see, the internal energy and the specific heat of the ensemble of oscillators 
do not depend on the dissipative media, at first order in the dissipation parame- 
ter. We must remain here that the whole system is made up of the crystal and 
the medium. So, in the above quantities one needs to add the contribution coming 
purely from the medium. On the other hand, we have used the classical canonical 
statistical partition function and the Debye's model for our study of the thermo- 
dynamic quantities of the system. However, this is somewhat a little bite incorrect 
since Debye's model works fine to relative low temperatures, and classical canonical 
partition function is expected to work fine to relatively high temperatures. Thus, 
let us make everything consistent by using quantum canonical partition function. 
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5. Quantum approach 

In this case, using (7), (20a), the condition (20b), and the same hypothesis of con- 
tinuity in the frequencies, one has the following expression 

3N / \ 
\nZ = haft E e_/3£nK) 

j=l \ n J 

= f> (]T e- 13 ^^) « [°° In (j2 e-P E ^) g{u) dcu 
j=1 



n / " u \ n 



Now, using (13) at first order in the dissipation parameter, one has 

1 - X 2 <r&*> ' 



(23) 



S.-WW (24o) 



where X a has been defined as 

Thus, Eq. (23) is written in the following way 



hiZ = -=-/ w 2 ln — ^To — AS- dw • (25) 



Therefore, the internal energy, the entropy, the specific heat, the equation of state 
and the free energy of the system are given by 

U = ^-L- [ uj 2 (Hlu — fiw a 7r) coth ( _ P^a^ X ^ 



9iV£; 



ln 1 - Xle-^ + 2 {TlU) ~ ^ C ° th H 2—, 



cv 



9iV 



j d uj 2 (huj — fiuJair) duj 



AkT^% Jo sinh 2 - ^) ' 



(27) 
(28) 



p - —In | 1 / u, 2 ln ^ A « e ^ 



l-A^e-/ 3 ^/ vu% Jo \1-X 2 a e-^ 



du: , (29) 



and 



9iV£;T , 

' oj ln 



w D Jo 



1 - A 2 , 



doj . 



(30) 



8 



Fig. 2 shows the variation of Cy/Nk as a function of the temperature for several 
values the dissipative parameter w Q . The effect at first order of the dissipation is to 
increases the specific heat a low temperatures. 



5. Conclusions 

We have used the phase space (x, v) to study the statistical properties of an en- 
semble of harmonic oscillators within a dissipative medium, where its effect on the 
oscillators is to create a linear velocity depending force. We have made the study at 
first order in the dissipation parameter for classical partition function and quantum 
partition function, taking the Debye's model of solids as an example for possible 
applications. The classical canonical partition function lead us to have an internal 
energy and specific heat which are independent on the dissipation. However, quan- 
tum canonical partition function brings about dependence on the dissipation for all 
thermodynamical variables of the system. It is our feeling that the use of the space 
(x, v) for statistical physics studies has less restrictions that the space (x, p). 
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Figure Captions 



Fig. 1 Function I(u a i pha /uj). 

Fig. 2 Specific heat with 9d = hojo/k = 150 K (corresponding to Sodium). 
(1): 150u a Tr/u D = 0.0, (2): 150u; Q 7r/u;o = 0.5, (3): 150w a 7r/w D = 1.0, (4): 
150w a 7r/u; D = 2.0, (5): 150w a 7r/w D = 5.0, (6): lbQuj a Ti/oj D = 10.0 
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